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Abstract 

The classical ferromagnet is analyzed in two dimensions via a mean-field treatment of 
the CP 1 model valid in the limit of low temperature and of weak magnetic field. At 
high skyrmion densities, we find that the magnetization begins to decrease linearly with 
temperature, but it then crosses over to a Curie-law at higher temperature. The comparison 
of these results with recent Knight-shift experiments on the quantum Hall state is fair. 
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It is well known that the continuum ferromagnet supports topologically non-trivial 
spin textures called skyrmions in two dimension. 1-3 In physical terms such skyrmions rep- 
resent domain-wall loops within which the magnetization is reversed. The scale invariance 
characteristic of the Heisenberg model in two dimensions indicates that the correspond- 
ing skyrmion solutions have no preferred size, nor do they interact, at zero temperature. 
Quantum fluctuations in the form of dynamical spin- wave excitations, however, break this 
scale invariance, which means that skyrmions interact at non-zero temperature. Similar 
effects arise in the case of the two-dimensional (2D) quantum antiferromagnet, but at zero 
temperature as well. 4 

Recent theoretical and experimental work indicates that skyrmion spin textures are 
tied to excess electronic charge with respect to a filled Landau level in the quantum Hall 
effect, 5-7 at which point the groundstate is a 2D ferromagnet. Motivated by this as well 
as by general theoretical considerations, we shall study here the magnetic properties of a 
gas of skyrmions with a net topological charge density. Like it's electronic counterpart, 
the total skyrmion number is necessarily quantized in integer values due to the topological 
nature of the 2D ferromagnet. 1-3 It is also a conserved quantity, which justifies the restric- 
tion to constant skyrmion number assumed throughout. The magnetic correlations, as well 
as the magnetization, are calculated as a function of temperature and of external magnetic 
field via a meanfield analysis of the CP 1 model for the classical ferromagnet. 2 ' 3 ' 8 Quantum 
effects 9 related to dynamical spin-wave excitations are therefore neglected. Notably, we 
obtain the following results in the limit of high skyrmion density: (a) the magnetic corre- 
lation length is set by the inter-skyrmion separation; (b) the paramagnetic susceptibility 
displays a Curie law at all temperatures due to the disordering effect of the skyrmion gas; 
and (c) the magnetization in fixed external field crosses over from a low-temperature linear 
decrease to the above Curie-law tail at high temperature. Although the latter cross-over 
is weak (see Fig. 1), the agreement between the theoretically predicted cross-over temper- 
ature with that observed in recent Knight-shift measurements 7 on the quantum Hall effect 
is surprisingly good. 

CP 1 Model. We now consider the thermodynamics of the classical 2D ferromagnet 
in the presence of a net skyrmion density. The energy functional of the corresponding 
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continuum Heisenberg model may be written in the CP 1 form as 

E = 2s 2 J J d 2 r\{V -ia)z\ 2 , (1) 

where z = (z-, z+) is a complex valued doublet field constrained to unit modulus, 

zz = 1, (2) 

and where the vector potential a is tied to the z fields by 

a = i(Vz)z. (3) 

Here s denotes the microscopic spin. The z fields are related to the normalized magneti- 
zation, m, of the Heisenberg model by 

m = zaz, (4) 

where a denotes the Pauli matrices, while the spin stiffness (1) is given by p s = s 2 J. 
Using Eqs. (2)-(4), it can be shown that the constrained (|m| 2 = 1) energy functional E = 
\p s J d 2 r\Vm\ 2 is equal to that of the CP 1 model (l). 2 ' 3 ' 8 Similarly, it can be shown that 
the skyrmion density is equal to the density of fictitious magnetic flux, (d x a y — d y a x )/2ir. 

To proceed further, we now write the Gibbs distribution, Z = J T>m8{\m\ 2 — l)e -E / fcBT ,| 
in the presence of external magnetic field H along the z-axis as 

Z = j VzVzVaVXexp^-2(3 j d 2 r{\ (V - ia)z\ 2 + i\(\z\ 2 - 1) + b z {\z-\ 2 - \z + \ 2 )} j, (5) 

where we define [5 = s 2 J/k B T and b z = nE z /4s 2 J, with spin density n and Zeeman energy 
splitting per spin E z = gpsH. Here integration over the Langrange multiplier field A(r) 
enforces constraint (2), while integration over the now unconstrained vector potential field 
a in the Coulomb gauge, V • a = 0, recovers constraint (3). We now make the basic 
approximations of the paper, which are (i) to enforce constraint (2) only on average over 
the entire 2D plane, and (ii) to neglect spatial fluctuations in the skyrmion density; 3 
i.e., the z fields above are integrated out first with the presumption of a homogeneous 
Langrange multiplier field A(r) = Ao and of a homogeneous fictitious magnetic field bs = 
(d x a y — d y a x ). The latter is of course equal to the net skyrmion density multiplied by 



a factor of 2n. The saddle-point condition -^InZ = is equivalent to the meanfield 
constraint 1 = (z-Z-) + (z + z + ), where 

(z±z±) = (2/3)" 1 y- 1 tr[-(V - id) 2 + i\ =F bz}- 1 (6) 

are the respective averages. But since the spectrum of the Hermitian operator —(V — id) 2 
is just that of Landau levels with energies 2bs(n+ \), each with a degeneracy per area V 
of bs/2ir, we obtain the meanfield equation 

2^/3= g 2hsn + i -2 + hz + 2 fl obsn + t-z-bz (?) 

for the correlation length £ set by £ -2 = zAo + 65. Similarly, we obtain the meanfield 
expression 
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for the normalized magnetization m z = (z + z + ) — (z-Z-) along the z direction. 

Consider now the ferromagnetic regime, 65 — > 0, where the concentration of skyrmions 
is dilute with respect to the magnetic correlation length, but is not zero. The sums in Eqs. 
(7) and (8) may then be converted into integrals, resulting in 2n/3 = In [/co(£~ 4 — ^z) -1 ^ 4 ] 
and m z = (27r/3) _1 ln [(£ -2 + bz)/(i~ 2 — bz)] 1 / 4 respectively. Here ko represents the 
momentum cut-off corresponding to the former integral. We therefore obtain the fa- 
miliar result £ = (£^~ 4 + frfz) -1 / 4 for the correlation length, 3 where the zero-field fer- 
romagnetic correlation length £ = k ~ 1 e 2nf3 diverges exponentially at low temperature. 
The magnetization is paramagnetic in the limit of weak external field, bz — > 0, fol- 
lowing m z = Xzbz with \z = (47r/9) _1 £ 2 . In general, it is given by the expression 
m z = 1 — (27r/3) _1 ln[/co(£ -2 + &z) -1 ^ 2 ], which means that the normalized magnetiza- 
tion begins to saturate logarithmically as external field increases. The continuum limit 
then breaks down at kobg 1 ~ 1, where m z is pathological and exceeds unity. 10 The 
present results correspond to those of the quantum ferromagnet in the renormalized clas- 
sical regime. 9 The former high-field catastrophe is avoided in that case, however, where 
the system evolves into a quantum activated regime instead. Also, it can be easily shown 
that (z±(f)z±(f')) ~ exp[— \f — r'\{^~ 2 =f bz) 1 / 2 ] at long distance \r — f'\, which by Eq. 



(4) indicates that (£ -2 — 6z) -1 ^ 2 gives the magnetic correlation length. Last, if 5bs{r)/2it 
denotes local fluctuations in the skyrmion density, while 5X(r) denotes the correspond- 
ing fluctuations of the Langrange multiplier field, then it can be shown that the Gibbs 
distribution associated with such fluctuations is given by 3 ' 8 

Z 2 = J D&^Aexpj-i J d 2 r[x(Sbs) 2 + 2ia(Sbs)(SX) + e(VSX) 2 ]^ (9) 

in the long- wavelength limit (relative to the inter-skyrmion separation), where x ~ £o> 
a = 2(3 /bs and e = (3/b 2 s (see the Appendix). Hence, the direction of steepest descent 
for gradients of the Lagrange multiplier field is along the real axis in the presence of a 
net skyrmion density. Notice also that homogeneous fluctuations (VSX = 0) of the Lan- 
grange multiplier field are in general only marginally stable in such case. Fluctuations, 
V~ 1 ((J d 2 r 5bs/2n) 2 ) ~ £q 2 , in the total skyrmion number are exponentially suppressed 
at low temperature, on the other hand. The latter thermally activated temperature de- 
pendence is then of course entirely consistent with the fact that Es = 4%s 2 J is the energy 
cost of a single skyrmion. 1 

Consider next the skyrmion-rich limit, bs — > oo, where the separation between neigh- 
boring skyrmions is much less than the ferromagnetic correlation length, £ - I* 1 that case 
the n = terms corresponding to the lowest Landau level dominate the sums in Eqs. (7) 
and (8), which results in the meanfield equations 2n(3 = &s£ 2 /(l — b 2 z ^ 4 ) and m z = ^ 2 bz 
respectively. Solving the former relationship explicitly, we obtain the simple expression 

mz = (l + t 2 ) 1 /2_ t - (1Q) 

for the normalized magnetization, where 

t = ^L^ (ii) 

E z nn v ' 

is the reduced temperature. The magnetization therefore decreases linearly from unity at 
low temperature and/or high magnetic field, while it follows a Curie law, m z = at 
high temperature and/or low magnetic field. (It is interesting to note that a similar H/T 
dependence is expected for the magnetization of the quantum ferromagnet in the high- 
temperature quantum-critical regime, yet in the absence of a net skyrmion density. 9 ) The 
effective area As(T) of a single skyrmion can be obtained from the comparison of Eq. (10) 



at low temperature with the naive formula m z = 1 — 2As(&s/27r) for the magnetization 

obtained by counting the number of reversed spins, S = nAs, per skyrmion. This yields an 

effective area of As = {ksE I Ez)n~ x that results from the opposition of the Zeeman energy, 

which tends to decrease the number of reversed spins, with the entropic part of the free 

energy, which tends to increase the number of reversed spins. The cross-over between these 

two regions naturally occurs at t = 1 (see Fig. 1), where the effective size of the skyrmion is 

by definition on the order of the inter-skyrmion separation. Also, the correlation function 

for the z fields has the gaussian form \(z(r)z(f '))| 2 = exp( — \\f — f '\ 2 bs) in above lowest- 

Landau-level approximation, which means that the magnetic correlation length is limited 

by the inter-skyrmion separation for all values of external magnetic field, bz, in this regime. 

This magnetic correlation length therefore also limits the size of a single skyrmion to the 

inter-skyrmion separation. Last, the long- wavelength fluctuation corrections to the present 

saddle-point (9) are identical to that of the skyrmion-poor regime, with the exception that 

X = 2f3/bs = o (see the Appendix). This means that fluctuations, V A_1 ((J d 2 r 5bs /27r) 2 ) ~ 

X -1 , in the total skyrmion number of the system are much larger in this regime. 

Quantum Hall Effect. Let us apply the above results to the quantum Hall effect in 

the vicinity of unit filling, where it has been suggested that the skyrmion density bs/2n 

is equal to the excess charge density 5n; 5 ~ 7 i.e., the electronic filling factor, v = 2-kIqTi, is 

related to the skyrmion density by \v — 1| = l$bs, where Iq denotes the magnetic length in 

fixed magnetic field. Using these relations, we then obtain the simple expression 

\v — 11 kuT , . 

t = 2 l -—^^r- (12) 
v hj Z 

for the reduced temperature (11). This immediately defines a cross-over temperature 
between the low-temperature regime with linearly decreasing magnetization and the high- 
temperature Curie-law regime (see Fig. 1) given by ksT* = ||1 — is~ 1 \~ 1 Ez, which is 
notably independent of the spin stiffness. Recent Knight shift measurements exhibit a 
temperature dependence that is qualitatively similar to that predicted by Eq. (10), yet 
with a more pronounced cross-over. 7 Such measurements were conducted in magnetic fields 
of H = 7T with a Zeeman energy of Ez — 0.2 meV, and at filling fractions of v = 0.88 
and v = 1.2. The present theory for the skyrmion rich regime then predicts cross-over 
temperatures of T* = 8K and T* = 7K, respectively. The sharp cross-overs observed 
experimentally at roughly 9 K in both cases compare quite well with these estimates. But 
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do such Knight-shift measurements actually lie within the skyrmion-rich limit discussed 
here? Inspection of Eq. (8) for the magnetization indicates that the skyrmion density 
must satisfy bs > max(bz, £,q 2 ) for this to be the case. Presuming that the momentum 
cut-off is set by the magnetic length, 11 ko ~ Iq , then the latter is equivalent to satisfying 
conditions (a) bs > bz, and (b) |1 — v\ > e ~ E s/k B T ^ w j iere e s = \-j\p s is the energy cost 
of creating a skyrmion in the ferromagnet. Standard estimates for the spin-stiffness p s 
of the quantum Hall state at unit filling 5 ' 9 place this energy at Es/ks ~ 40 K in a 7T 
field, which means that condition (b) is easily satisfied at temperatures much less than 
this scale. Also, employing Eq. (11) at the cross-over (t = 1) implies that condition (a) 
is satisfied as well, since T-m ~ 10 K falls substantially below the previous scale. Last, the 
present classical treatment of the ferromagnet is generally valid in the renormalized classical 
regime, 9 Ez < ksT < E$, which spans the majority of the experimental temperature 
range. 7 

In conclusion, although the magnetization (10) obtained here for the classical 2D 
ferromagnet in the presence of a net skyrmion concentration compares favorably with re- 
cent Knight shift experiments on the quantum Hall state, it consistently overestimates the 
latter. 7 From the theoretical side, a number of effects could be responsible for this discrep- 
ancy. First, it has been noted that fluctuations with respect to the present meanfield treat- 
ment of the classical ferromagnet are relatively large in the presence of skyrmions. More 
generally, it is also known that fluctuations in the skyrmion density act to confine the z and 
z fields in strictly two dimensions. 8 Such effects have not been accounted for here. Second, 
the inclusion of quantum mechanical spin-wave excitations will further decrease the mag- 
netization. Last, both quantum fluctuations 4 as well as the direct Coulomb interaction 12 
tend to crystallize the skyrmion gas. Accounting for all of the above-mentioned effects 
remains as a problem for the future. 

The author has benefited from numerous discussions with L. Brey and A. Somoza. 
This work was supported in part by the Spanish Ministry for Science and Education, as 
well as by National Science Foundation grant DMR-9322427. 
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Appendix 



We shall now compute the fluctuation corrections (9) to the present meanfield approx- 
imation. If we take the Landau gauge, a = (0, bsx), then the eigenstates of the Hermitian 
operator —(V — ia) 2 + iXo have the usual form (r\n,q) = (x'\n)e~' iqy , with eigenvalues 
e n = 2bsn + £ -2 . Here \n) denotes the eigenstates of the corresponding Harmonic os- 
cillator while x' = x — bg q. If we now identify the variation in the vector potential by 
Sbs — d x 5a y — d y 5a x , then fluctuations correct the Gibbs distribution by a factor of 

exp \-\v- 1 £{ PS* + n,i(*)]M*)M-*) + 2<[nj(fc) + u-(k)]d aj (k)d\(-k)+ 
V k 

+ [U+ (k) + U^ (k)]SX(k)SX(-k)}^ (Al) 

where Ilfj(k), U.Q Q (k) and are respectively the high-temperature limits of the Kubo 
formula, the Lindhard function and the hybrid Kubo-Linhard function corresponding to 
the z± field. Also, 5di(k) and 8X(k) denote the appropriate Fourier components. One may 

— * 

presume that k be parallel to the x-axis without any loss of of generality, in which case 
the polarizabilities are given explicitly by 

n±/^_TMV^ 25 ij t/-i V" ( n > K, q) (n', q\vj(-k)\n, g) 

n,q n^n' ,q 



n ± f%\ = y-i (n,q\vj(k)\n',q) (n\q\e ikx \n,q) 



£ n =F bz e n > =F b z 

n=£n,q 



(A3) 



n± \y-i {n,q\e lkx \n',q} (n',q\e lkx \n,q) 

00 • ' ' f-f s„ : by, e n > | b x 

n=£n,q 

where v(k) = 2e l ^ kx (iV + a)e l ^ kx is the modified velocity operator. In general, these 
polarizabilities have the form n^ (/c) = e±k 2 , U^y(k) = a±ik and n^(fc) = x±k 2 in 
the long- wavelength limit, while the remaining ones are identically zero due to gauge 
invariance. A direct application of (A3) and (A4) then yields a± = b^V -1 ^EnT^) -1 
and e± = b~^ 2 \V~ l J2 n q ( e n Tbz) -1 - Employing mean-field Eq. (7), we then obtain 
that the net fluctuation susceptibilities (9) are given by a = a + + cr_ = 2/3/bs and 
e = e + + e_ = f3/b 2 s in general. Similar calculations for (A2) also demonstrate ultimately 



that x = 2/3/6s = a in the skyrmion-rich limit, bs — > oo. This identity is suggested by the 
present mean-field analysis, which depends crucially on the combination £ -2 = iAo + bs 
[see Eq. (9)]. 
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Figure Caption 

Fig. 1. Shown is mean-field result (10) for the magnetization vs. temperature of the classical 
2D ferromagnet in the limit of high skyrmion density. The pronounced Curie-law tail 
notably weakens the cross-over at T*, which should lie well below 4itp s /kB- 
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